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Abstract 

We investigate path integral quantization of two versions of unimodular grav¬ 
ity. First a fully diffeomorphism-invariant theory is analyzed, which does not 
include a unimodular condition on the metric, while still being equivalent to 
other unimodular gravity theories at the classical level. The path integral has 
the same form as in general relativity (GR), except that the cosmological con¬ 
stant is an unspecified value of a variable, and it thus is unrelated to any coupling 
constant. When the state of the universe is a superposition of vacuum states, the 
path integral is extended to include an integral over the cosmological constant. 
Second, we analyze the standard unimodular theory of gravity, where the metric 
determinant is fixed by a constraint. Its path integral differs from the one of 
GR in two ways: the metric of spacetime satisfies the unimodular condition only 
in average over space, and both the Hamiltonian constraint and the associated 
gauge condition have zero average over space. Finally, the canonical relation 
between the given unimodular theories of gravity is established. 


1 Introduction 

The idea of unimodular gravity is nearly as old as general relativity (GR) itself. Orig¬ 
inally, Einstein considered the unimodular condition [l|, 

( 1 - 1 ) 

as a convenient way to partially fix a coordinate system in GR, which simplihes the 
calculations in certain situations. Later on, unimodular gravity has also been consid¬ 
ered as an alternative theory of gravity closely related to GR, which was first suggested 
in . The definition of unimodular gravity is usually based on the invariance under a 
restricted group of diffeomorphisms that leave the determinant of the metric invariant, 
so that the determinant of the metric can be set equal to a fixed scalar density Cq, 

= ^05 ( 1 - 2 ) 

which provides a fixed volume element in spacetime. We consider a theory based on 
the condition fll.2j) . and on the associated restricted group of diffeomorphisms, as the 
hrst example of unimodular gravity. 
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Fully diffeomorphism-invariant extensions of unimodular gravity exist as well, which 
also involve a condition on the determinant of the metric such that the right-hand side 
of the condition fll.2p is replaced with a scalar density held. The most prominent the¬ 
ory of this kind is the Henneaux-Teitelboim theory , where the unimodular condition 
sets ^/^ equal to the divergence of a vector density held. 

It is well known that classically unimodular gravity produces the same physics 
as GR with a cosmological constant. The held equation for the metric is either the 
traceless Einstein equation or, thanks to the Bianchi identity, the Einstein equation 
with a cosmological constant |^. The diherence is that the cosmological constant of 
unimodular gravity is a constant of integration, rather than a coupling constant. Since 
the value of the cosmological constant is unspecihed and unrelated to any coupling 
consta nt, p roblems associated with the cosmological constant have been reconsidered 


see 


j^lO for reviews) 


Quantum corrections to the energy-momentum tensor of matter, which are 
of the form where G is a constant over spacetime, do not contribute to the 

traceless held equation for the metric in unimodular gravity. In particular, vacuum 
huctuations in the trace of the energy-momentum tensor of matter do not ahect the 
metric. This well-known feature of unimodular gravity has been recently revisited 
via an explicit calculation of one-loop corrections llj. Since a small nonvanishing 


cosmological constant is required, the full Einstein equation and an associated action 
need to be considered. There the vacuum corrections are absorbed into the arbitrary 
cosmological constant, whose value should be specihed experimentally. However, this 
does not solve the cosmological constant problem. Unimodular gravity faces a similar 
OToblem with the renormalization or hne tuning of the cosmological constant as GR 
B0SE3- The expression for the vacuum energy generated by the quantum huc¬ 
tuations is highly dependent on the details of the ehective description, in particular 
on the chosen Wilsonian cut-oh scale B]. Therefore we do not consider the vacuum 
energy problem in this paper. Instead we concentrate on the formal diherences be¬ 
tween the unimodular gravity theories and GR at the quantum level. In other words, 
our treatment assumes that somehow the observed cosmological constant A will be 
stabilized against vacuum corrections. 

Predicting or deriving the observed value of the cosmological constant is a hard 
problem as well. A highly speculative but interesting attempt to address this problem 
in unimodular gravity has been made in [HQ, where an integral over the cosmo¬ 
logical constant was included into the path integral. We will show how this argument 
can be derived in a straightforward way, when a new action for unimodular gravity is 
introduced. Problems associated with the given argument are also discussed. 

Gonventionally, the idea of unimodular gravity has been to impose a condition 
on the determinant of the metric, e.g., fll.2p . In comparison with GR, making the 
cosmological constant an arbitrary constant of integration can be regarded as the key 
feature of unimodular gravity. In order to achieve it, however, there is no need to 
constrain the determinant of the metric. We consider a fully diffeomorphism-invariant 
theory (see fl2.13ll h which has recently appeared in the context of gravity with mimetic 
dark matter [1^, where an additional scalar held was also included to describe the 
mimetic matter. The given theory is no longer unimodular in the sense that there is 
no condition on the determinant of the metric, but we will establish how the theory 
is canonically related to the conventional unimodular theories of gravity. 

It has been argued that unimodular gravity can offer a new perspective on the 
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problem of time in quantum gravity and cosmology [4, Il3|, |ly, |17[. Since the bulk 


part of the Hamiltonian of unimodular gravity is nonvanishing, and the four-volume 
provides a cosmological time, an analogy of the Schrodinger equation exists, and hence 
quantum states of the universe can evolve in terms of a global time. On the other 
hand, it has been concluded that unimodular gravity cannot solve the problem of 
time in quantum gravity [l^, since the four-volume labels only equivalence classes of 
hypersurfaces separated by a zero four-volumes. 

Since all versions of unimodular gravity must be classically equivalent to GR, 
quantization of each version of unimodular gravity can be regarded as a potential 
quantization of GR. Hence it is necessary to understand how the different versions of 
unimodular gravity differ from each other and from GR at the quantum level. The 
equivalence of GR and unimodular gravity was recently discussed in 1^, concluding 
that the equivalence can be retained at quantum level when the UV extension of uni¬ 
modular gravity is performed appropriately. We will see that the form of the path 
integral depends on which version of unimodular gravity is chosen. 

Path integral quantization of the Henneaux-Teitelboim version of unimodular grav¬ 
ity has been considered previously in Q (see also |^), where the unimodular con¬ 
dition was shown to be imposed locally in the quantum theory. In this paper, we 
study the path integral quantization for the two other versions of unimodular gravity 
discussed above. The results are compared to both GR and the Henneaux-Teitelboim 
theory. In the fully diffeomorphism-invariant theory (see fl2.13jl for action), the path 
integral has the same form as the one of GR with a cosmological constant, but the 
value of A is an unspecihed constant value of a variable. Two approaches regarding 
the interpretation of the cosmological constant are considered: either (i) the effective 
value of the cosmological constant is hxed by the physical boundary conditions of the 
path integral, or (ii) the state of the universe is taken as a superposition of states with 
different values of A, and consequently the path integral includes an integral over A. 
In the latter a ppr oach, we derive the path integral in the form originally proposed in 
(l^ (see also [l3|, [l^). In the theory with a hxed metric determinant (see fl2.6p for 
action), the unimodular condition fll.2p is found to be imposed in average over space, 
but not locally. 

In Sect. [2] we present the different actions of unimodular gravity which are relevant 
for this paper, and discuss how the (classical) actions are related to each other. Sec- 
tionOis devoted to the canonical path integral quantization of the fully diffeomorphism- 
invariant theory of unimodular gravity. In Sect. Hlthe same is achieved for the con¬ 
ventional version of unimodular gravity with a hxed metric determinant. Section |5] 
establishes the canonical relation of the theories. The results are discussed in Sect. El 


2 Three versions of unimodular gravity 


2.1 Unimodular gravity with a fixed metric determinant 

Gonventionally, the held equations of unimodular gravity are obtained from the Einstein- 
Hilbert action under a restricted variation of the metric that preserves the deter¬ 
minant of the metric, 

^ -V^ = 0, ( 2 . 1 ) 




flU 
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where g = det g^iy. Since the metric transforms under an infinitesimal diffeomorphism, 




( 2 . 2 ) 


as 

fiiv + vitii (2.3) 

the unimodular condition fl 2 .ip requires that 

= 0, (2.4) 

i.e., 

= 0. (2.5) 

These transformations are often referred to as transverse diffeomorphisms or volume¬ 
preserving diffeomorphisms. However, the name transverse diffeomorphisms (TDiff) 
is sometimes reserved for the transformations that satisfy the noncovariant condition 
= 0 . In order to avoid any confusion, we shall refer to the given transforma¬ 

tions fl2.ip - fl2.5p as metric determinant-preserving diffeomorphisms. 

One way to define unimodular gravity is to introduce the unimodular condition 
fll.2p into Einstein-Hilbert action as a constraint multiplied by a Lagrange multiplier 

A, 


S\jG[gfiu,K'^] = J - A (\/^- 

+ - (f d^X^/\^\K: + Sra[g^^u, 

^ JdM 


( 2 . 6 ) 


where eo is a fixed scalar density, such that eod^x defines a proper volume element, the 
gravitational coupling constant is denoted as k = IbvrG, and is the action for the 
matter fields (denoted collectively by \h) which are coupled to the metric in the same 
way as in GR. In the surface integral over the boundary dAi of spacetime, 7 is the 
determinant of the induced metric on dAi, and K, is the trace of the extrinsic curvature 
of the boundary. The boundary term is included as in GR, so that the variational 
principle for the action is well dehned without imposing boundary conditions on the 
derivatives of the metric]^ The full diffeomorphism invariance of GR is lost due to the 
presence of the hxed volume element eod^x. The action fl2.6p is invariant under the 
metric determinant-preserving diffeomorphisms. We shall refer to the theory defined 
by fl 2 . 6 p simply as unimodular gravity (UG). 

An unrestricted variation of g'^^ gives the Einstein equation 


Rfiv 


—^gfj.u = ^2/XI/. 


(2.7) 


The variation of A gives the unimodular condition fll.2p . The field equations for matter 
fields are identical to those of GR. The unimodular condition fll.2p ensures that fl2.ip 
holds. The energy-momentum tensor of matter is defined in fl2.7l) as usual, = 
We assume that the action for matter is diffeomorphism-invariant, so that 

“^When we write about boundary conditions without specifying their nature in the canonical for¬ 
malism, we refer to both the initial conditions and the conditions on the spatial boundary. Likewise 
in Lagrangian formalism we refer to conditions on the boundary of spacetime. 


4 









energy-momentum is conserved, = 0. Then we take the divergence of fl2.7p and 

obtain 

V^A = ^ = 0, (2.8) 

where the (contracted) Bianchi identity is used. Thus we see that A is hxed as a 
constant of integration, which we denote as A = -^A, where A is the cosmological 
constant. Inserting this into the held equation fl2.7p gives 

1 Ki 

Rfiu '^Rq^j.u + A^f^jy = (2-9) 

Compared to GR the only difference is that we are restricted to use coordinate sys¬ 
tems that satisfy fll.2p . In GR, the condition fll.2p can always be satished locally by 
choosing the inertial coordinates. Then every coordinate system obtained via metric 
determinant-preserving diffeomorphisms satishes fll.2p as well. 


2.2 Fully diffeomorphism-invariant unimodular gravity 

Extensions of unimodular gravity with full diffeomorphism invariance have been pro¬ 
posed as well. The nmst prominent theory is dehned by the Henneaux-Teitelboim 
(HT) action 


see 


22 | regarding the boundary surface term). 


A,r'^, T] = j d^x (v^ “ 

-F (f d^x (-\/yi\K: - Ar^r^^ -F T], 

JdM \ ^ / 


( 2 . 10 ) 


where is a vector density and is the outward-pointing unit normal to the boundary 
dA4. The held equations consists of the Einstein equation fl2.7p . the equation for the 
cosmological constant variable, 

V^A = 0, (2.11) 

a (fully diheomorphism-invariant) unimodular condition. 


V^ = d, 




( 2 . 12 ) 


and standard held equations for matter. 

The HT action 02 .101) can indeed be derived from the UG action 02.61) via pa¬ 


rameterization of the spacetime coordinates [18|. Parametrization of coordinates in 


a mechanical system is a well known method for obtaining a reparameterization in¬ 
variant action (see for a review). Parametrization of held theories was introduced 
later (see for a description). We treat the coordinates of the action 02.6p as four 
independent scalar variables X“(a;) that depend on the actual coordinates x^. One can 
think of this as a transformation —)■ X“(x). The Einstein-Hilbert and matter parts 
of the action 02 .6p are invariant under such transformation, but the part with a hxed 
volume element is not invariant, since it transforms as J d'^xeoX —)■ / (i^a;eoA|c?^X“|, 
where \d^X°'\ is the Jacobian determinant of the transformation. When we identify 
a vector density as we obtain the HT action 

O2.10p . It is clear that the HT theory is classically equivalent to the UG theory 02.6p . 
However, diherences are expected to arise upon quantization. 
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We consider an alternative action that is fnlly diffeomorphism-invariant and retains 
the classical eqnivalence with the other nnimodnlar theories, in particnlar with fl2.6p 
and fl2.1UI) . The action has been stndied in the context of gravity with mimetic dark 
matter [l5|, where an additional scalar held was also inclnded to describe the mimetic 
matter. The action is written (withont the scalar held) as 

[ d^x^(--X-V^V,x) 

Jm J ('213) 

d^x\/\^\lC + T], 

where the variable is a vector held. We shall refer to this theory as the fnlly 
diheomorphism-invariant nnimodnlar gravity (DUG). The action 02.131) is argnably 
the most transparent dehnition of snch a theory. The action 02.13|) consists of the 
Einstein-Hilbert action with a variable cosmological constant A, and a constraint term 
for A. The vector held acts as a Lagrange mnltiplier that ensnres V^A is zero in 

every direction, and thns A is a constant. The held eqnations consists of the Einstein 

eqnation 02.71) for the metric, the eqnation 02.lip for the cosmological constant variable 
A, an eqnation for the anxiliary vector held 



= 1 , 


(2.14) 


and standard held eqnations for matter. The nnimodnlar condition on the metric de¬ 
terminant, 01.2p or 02.12p . has been replaced with the condition 02.14p on the vector 
held. The vector held does not contribnte to the Einstein eqnation dne to Eq. 02.111) . 
In Sect. [21 we will show how the vector held can be eliminated from the Hamilto¬ 
nian formnlation while the canonical representation of diheomorphism invariance is 
retained. 

It is obvions that the DUG action 02.131) is closely related to the HT action 02.101) . 
An integration by parts in the term J d'^xXd^r^^ of the HT action, followed by a 
replacement of the vector density variable with a vector held variable, = \/—gV^, 
gives the action 02.13p . Hence it is clear that these theories are eqnivalent classically. 
However, the path integral for the action 02.13p will be shown to diher from the HT 
case signihcantly dne to the diherent choice of variable. 

The held eqnations for both the HT and the DUG theories are invariant nnder the 
shift 

T/iu T/iu + A —)■ A -|- G, (2-15) 


where G is a constant. Hence qnantnm corrections to the trace of the energy-momentnm 
tensor are absorbed into the variable A, whose valne is an arbitrary constant. The vari¬ 
able A will be shown to remain constant at qnantnm level in Sect. [3l 

There exist more versions of nnimodnlar gravity in addition to the three theories 
discnssed above; see, for example, 14, [13 for other actions. In this paper we will 
concentrate on the three theories dehned by 02.61) . 02.10p and 02.13p . 


6 




































3 Quantization of the fully diffeomorphism-invariant 
unimodular gravity 

3.1 Arnowitt-Deser-Misner decomposition of the action 

Spacetime is assumed to admit a foliation to a union of nonintersecting spacelike 
hypersurfaces. The hypersurfaces are labelled by a scalar t that is constant across 
each hypersurface. The future-pointing unit normal to is denoted by n^. The so- 
called direction of time vector satishes = 1. The metric has the signature 

and hence = —1. Each hypersurface is described by the induced 

metric on Sj, 

and by the extrinsic curvature tensor 

Kfj^iy = ^-|- (3-2) 

where we have dehned the acceleration vector of Eulerian observers as 


= n''Vi 


n 


/i- 


(3.3) 


Then we introduce the Arnowitt-Deser-Misner (ADM) variables. The scalar t is 
taken as the time coordinate. The unit normal to is written as 


n 


0 


1 

iV’ 


n 


iV* 

Iv’ 


(3.4) 


where N is the lapse variable and V* is the shift vector on S^. Latin indices {i,j ,...) 
range from 1 to 3. Now the metric takes the form 


goo = + NiN\ goi = Ni, gij = hij, (3.5) 


where Ni = hijNK The extrinsic curvature is written as 


Kij = {dthij - DiNj - DjNi ), (3.6) 

where D is the covariant derivative that is compatible with the metric hij on and 
is the inverse metric, hijh^^ = 5A. The trace of extrinsic curvature is denoted by 
K = 

The action is decomposed as follows. The metric determinant is given by 

= NVh, (3.7) 


where h = det hij. The scalar curvature is written as 


R = + 2V^(n^iL - a^), 


where the De Witt metric is dehned as 


gijkl ^ 


(3.8) 


(3.9) 
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and is the (intrinsic) scalar cnrvatnre of S*. The last term in fl3.8p is a total 
derivative which contribntes a bonndary term into the action. 

The vector held is decomposed into components tangent and normal to as 



(3.10) 

where 


±V^ = K, = 

(3.11) 

and the projection operator onto S* is dehned as 



(3.12) 


The gravitational part of the action fl2.13p is written in ADM form as 


SnvG[N,N\hij,X,Vr,,V\'^]= / dt / d^xNVh 


'St 


K 




A + 14VttA - V^diX 


+ Sl3 + Sra[g^j.u,'^], (3.13) 


where we denote D* = and 

VttA = ^ {dtX - N%X) , (3.14) 

and the bonndary contribution Sb is given as in GR, 

Sb = — [ d^x^/—'y (JC + r^n^K — r^a ^), (3.15) 

« Jb 

where B is the timelike part of the bonndary dAd. The snrface B is foliated to a nnion 
of two-dimensional snrfaces i?*, which come from the intersection of St and B. When 
the hypersnrfaces B and St are orthogonal, the snrface term 03.151) can be written as 

3 

= l [ dt [ d^xNy/^^^^K, (3.16) 

J Jst 

where a is the determinant of the induced metric on Bt, and is the trace of the 
extrinsic curvature of Bt in St. 


3.2 Hamiltonian analysis 


Hamiltonian analysis of unimodul ar g ravity in its different forms has been considered 
in several papers [i, 0, [H, 17, 18, 20, 26, 27|. Since the action 02.131) differs from the 
previous theories by lacking a unimodular condition and involving the vector held, we 
present a detailed Hamiltonian analysis. 


3.2.1 Hamiltonian and constraints 

We shall obtain the Hamiltonian and the full set of constraints for the action 02.13p . 
Here we consider pure gravity, since the matter sector is identical to that of GR, and 
in the end we include matter into the path integral in Sect. 13.31 
























First we introduce the canonical momenta tt^v, tt*, px, Pi and pn conjugate 
to N, N\ hij, A, id* and Vn, respectively. Since the action fl3.13p is independent of 
the time derivatives of the variables N, N^, Vh and Id*, their canonically conjugated 
momenta are primary constraints: 


tin ~ 0, TliTH 0, Pn ^0, Pi^i 0. 


(3.17) 


In addition, the dehnition of the momentum conjugate to A implies the primary con¬ 
straint 

Cx=px- ^/hVn ~ 0. (3.18) 

The momentum conjugate to the metric is dehned as 

The Hamiltonian is obtained as 


H= (fx{N'HT + N"'Hi + VNnN + v]^7li + VxCx + VnPn+v"Pi)+Hi3^, (3.20) 

JT,t 

where the so-called super-Hamiltonian and supermomentum are dehned as 

Ut = +x/hX + x/hV%X (3.21) 

Vh K 

and 

Hi = -2hijDkTi^^ + d^Xpx, (3.22) 

respectively, where we introduced the inverse De Witt metric as 

dijki '^(^hikhji T (3.23) 

and VN,v}^,vx,v\Vn are unspecihed Lagrange multipliers for the primary constraints. 
Regarding the surface terms, the analysis follows the standard set by jUj . The surface 
term in the Hamiltonian 03.201) is obtained as 

Hb, = -- [ + 2 [ (fxN^hijrkTi^\ (3.24) 

^ JBt JBt 

which is the same expression as in GR. The physical Hamiltonian is dehned with 
respect to a chosen reference background as ifphys = H — Ffref, where the Hamiltonian 
of the background is denoted as Ffref- The total gravitational energy of the system is 
the value of the physical Hamiltonian. The surface term 03.24p is given in a generic 
form that produces the correct expression of total gravitational energy for diherent 
reference backgrounds j^. 

We must ensure that every constraint is preserved under time evolution that is 
generated by the Hamiltonian 03.20^ . The preservation of ttjv ~ 0 is ensured by the 
Hamiltonian constraint 

nr^O (3.25) 

and the preservation of tt, ^ 0 is ensured by the momentum constraint 


n. 


0 . 


(3.26) 
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We can extend the momentnm constraint fl3.22jl with a term that is proportional to 
the primary constraint Pn so that the momentnm constraint generates spatial diffeo- 
morphisms on for all the variables that are involved in the constraints]^ For that 
reason we redehne 


Ui = -2hijDkTi^^ + diXpx + diVnPn ^ 0. (3.27) 

It is nsefnl to dehne global (smeared) versions of these constraints for calcnlational 
pnrposes: 

= [ d^x^T, <h[xi = [ (3.28) 

where ^ and y* are fnnctions on The preservation of the constraint pi ~ 0, 

dtPi = {Pi, H] K. -Ny/hdiX 0, (3.29) 

is ensnred by introdncing a new constraint, 

Ci = diX ^ 0. (3.30) 

This constraint implies that A is a constant across S*. We dehne the smeared form of 
Ci as 

^[x"] = [ d^xx'diX. (3.31) 

JEt 

This constraint is inclnded into the Hamiltonian with a Lagrange mnltiplier as 
The preservation of the constraint Ca ~ 0, 

dtCx = {Cx,H} ^ {CxMN]} - Vhvr, + {CxA<]] ~ 0, (3.32) 

is ensnred by hxing the Lagrange mnltiplier Vn of the constraint pn as 

Vn = -N + (3-33) 

The preservation of the constraint Pri ~ 0, 

dtPn = {Pn, H} ^ y/hvx ~ 0, (3.34) 

is ensnred by hxing the Lagrange mnltiplier vx of the constraint Cx as 

vx = 0. (3.35) 

Since the constraint Ci is inclnded into the Hamiltonian with a Lagrange mnltiplier, 
we can simplify the system by redehning T-Lt withont the part that is proportional to 
Ci- Now the total Hamiltonian is written as 

H = [ d^x {NA + N'^'Hi + vnt^n + vAi + x^'Pi + v^Ci) + HiSt, (3.36) 

JEt 

®We do not need to include a generator for the variables and pi, since the F*-dependent term 
in the Hamiltonian constraint (13.211) is proportional to the constraint (13.301) found below. 


10 










where we have defined the constraints 


n'T = nT-Pr. + ~ 0, (3.37) 

Ht = + \/hA 0, 

Vh K 

C' = C,-d, ^ 0, (3.38) 

and VN,vlf,v\v\ are unspecified Lagrange multipliers. 

What remains to be established is the preservation of the secondary constraints Rt, 
Ri and Cj under time evolution. The constraints Rt-, Ri, Ci have vanishing Poisson 
bracket with pn. The constraints Rt, Ri, Ci satisfy the following algebra: 

{RT[i],RTVl\} = [ d^x{^dip - r]di^)R^{Rj - pxCj - djVnPn), 

JTt 

mx%nT[^]} = RT[x%^], 

{Rt[^],C[x']} = 0 , 

mx%CW]}=C[x^d,p^]. (3.39) 

The first three Poisson brackets in fl3.39p are the familiar relations found in GR. The 
constraint in the right-hand side of the first Poisson bracket is just the momentum 
constraint of GR, Ri—pxCi — diVnPn = —‘^^hijDk'n^^ ~ 0. The last two Poisson brackets 
tell that Ci is preserved in time and that it transforms as a vector density under the 
spatial diffeomorphisms generated by the momentum constraint. The constraints ttat, 
TTj, Pi have vanishing Poisson bracket with every constraint. Thus all the constraints 
are preserved under time evolution. 

We can now see that all the constraints ("H^, "Hi, vtjv, C') in the Hamiltonian 
fl3.36p are first class constraints. The Lagrange multipliers in the Hamiltonian fl3.36p 
remain unspecified, until they are determined as a part of the gauge fixing procedure. 
In addition, p^ ^ 0 and Ca ~ 0 are the second class constraints. 

In order to clarify the nature of the constraint Ci, it is useful to decompose the 
variables X,px as follows: 


X{t,x) = Xo{t) + X{t,x), 


P\{'t,x) = 


x/h g 


d^xx/h 


P\{t) +Px{t,x), 


(3.40) 


where the zero modes describe the time dependent averages of A and px over space, 
Xo{t) = . ,3 / 7 - / d^x\fhX{t,x), pl{t)=f d^xpx{t,x), (3.41) 

J ^ (Jj tX/ IV J J St 

and the barred components have vanishing average values over space. 


/ d^xx/hX{t,x) = 0, / d^xpx{t,x) = 0. 

'st ist 


(3.42) 
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If the space Sj is infinite, the dehnition of the zero modes fl3.4ip has to be specihed 
more precisely. For example, in the asymptotically flat case, the spatial integrals would 
be dehned up to a hnite radius r in the asymptotic region, and hnally the limit r —)■ oo 
would be taken. In the dehnition of Aq the two inhnite integrals cancel out, since 
the asymptotic value of A must be a constant, so that the average value Aq remains 
hnite. The momentum p\ can be dehned to have such an asymptotic behavior that 
the dehnition of its zero mode remains hnite. Other scalar helds or scalar densities 
can be decomposed in a similar way. The zero modes satisfy the canonical Poisson 
bracket 

{Ao,pU = 1, (3.43) 

while the average free components satisfy 

= (3.44) 

JEt ^ 

and the Poisson brackets between zero modes and average free components are zero 

{Ao,PA(a:)} = 0, {A(a;),p°} = 0. (3.45) 

When A is decomposed, the constraint fl3.30p can be replaced with a local constraint 


A 0, 


(3.46) 


since c?jA = c?jA = 0 implies that A is constant over space and the zero average condition 
fl3.42p requires that constant to be zero. The corresponding hrst class constraint fl3.38p 
is replaced with 


C = 



(3.47) 


where the overline denotes a component whose integral over space vanishes. The 
purpose of the decomposition fl3.40p of the cosmological variable is to separate the 
perturbative component A that vanishes due to the constraint fl3.3Up . The average 
component Aq is left unconstrained. 

The total Hamiltonian fl3.36p is rewritten as 


H= f d^x [NH't + N'^'Hi + vnt^n + + vxC + v^'pi) + (3.48) 

where the variable A in the Hamiltonian constraint fl3.37p is replaced with its zero 
mode Aq. Next we consider gauge hxing and simplihcation of the Hamiltonian via 
elimination of some variables. 


3.2.2 Gauge fixing and the second class constraints 

Each of the hrst class constraints generates a gauge transformation. The constraint Pi 
generates a gauge transformation of the vector W as 

= |^^ y fcVij = £\ (3.49) 

where is an inhnitesimal gauge parameter. This means that W can be hxed through¬ 
out spacetime as a gauge choice. We choose the gauge hxing condition as = 0. We 


12 




















can further simplify the system by considering the gauge symmetry that is associated 
with the constraint fl3.47p . The constraint fl3.47p generates the transformation of the 
average free momentum px, 


Spx 



-e, 


(3.50) 


where the inhnitesimal gauge parameter £ is a now a scalar density whose integral over 
St vanishes. Equivalently, the constraint fl3.38p generates the transformation 


Spx = 



die" 


(3.51) 


where the integral of the component of the inhnitesimal gauge parameter e" in the 
direction of the outward-pointing unit normal r* to the boundary of St is zero, so 
that dP’xdie" = dPxrie" = 0. The corresponding gauge freedom can be hxed by 
setting Pa = 0 . 

Now we have the set of second class constraints {Cx,Pn,V",Pi,\,px)- The second 
class constraints can be set to zero strongly, if we replace the Poisson bracket with the 
Dirac bracket. In this case the Dirac bracket is equal to the Poisson bracket. Then we 
can eliminate six canonical variables iVn,Pn,V",Pi,X,px) by using the constraints 


A 


Pa 

eh’ 


0 . 


0 , Pi = 0 , 

0 , Px = 0. 


The Hamiltonian fl3.48p is written as 


(3.52) 


H= d^x [NKt + N^Tii + vnttn + v"xf7ii) + HiSt, (3.53) 

where 

Kt = + VhXo ^ 0, (3.54) 

x/h le 

and 

ni = -2hijDk7r^'" ^ 0. (3.55) 

This is the Hamiltonian of GR with a time dependent cosmological constant Aq. How¬ 
ever, it is evident that Aq is a constant in time as well, since the Hamiltonian is 
independent of 

dfXo = {Xq, H} = 0. (3.56) 


The value of Aq is set as a part of the initial value data on the initial Cauchy surface, 
say So at f = 0. The momentum p° evolves monotonically, dtp^ = —N\/h, and it is 
not involved in the actual dynamics of the system. The physical degrees of freedom 
consist of the two standard modes of gravity, plus the nondynamical zero mode that 
provides the cosmological constant. 

Gauge hxing the diffeomorphism invariance is done exactly as in GR by introducing 
appropriate gauge conditions for the generators Ht, 'Hi, tin, tt*. 
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3.3 Path integral 

3.3.1 Canonical path integral and possible ganges 

The canonical Hamiltonian for the gravitational part of the action fl2.13jl is written as 


He = / d^x [NUt + N^Hi) + Hb„ (3.57) 

where the Hamiltonian and momentnm constraints are defined in fl3.21l) and fl3.22l) . 
and the bonndary term in 03.241) . The second class constraints are 

Pn ~ 0, Cx=Px- VhVn ^ 0. (3.58) 

The hrst class constraints are 7rAr!^0,7ri~0,pi~0,'Hi~0 and 

= 'Ht - P-n + -^j^Vr^Prv ~ 0 , 

We introdnee gange hxing conditions as 

= N - f ^0, a* = W - ^ 0, W 0, 

TT*^) ^0, Pa ~ 0, (3.60) 

where / and /* are fixed fnnetions (or constants), snch that / > 0, while the conditions 

can depend on both hij and presnming that x^ depends on tt®-^ linearly or not 

at all. The fonr gange conditions x^ have to be independent, so that they £x fonr 
components of the variables hij,p'^^. Fnrthermore, it is convenient to reqnire that 

{x^x"} = 0. (3.61) 

The generator C' exhibits a nonlocal linear dependence over the spatial hypersnr- 
face, since the spatial integral of the generator vanishes by definition. The correspond¬ 
ing gange condition [px ~ 0) has a similar nonlocal linear dependence. Qnantization of 
gange theories with linearly dependent generators is disenssed in(see also (i^b 
There we show that the following path integral is obtained when certain additional 
gange conditions are imposed on the ghost fields associated with the generator C'. 

Since the first class constraints have vanishing Poisson brackets with every con¬ 
straint except the gange conditions, we nse the Faddeev formnla for the fnnctional 
determinant of constraints in the path integral. Furthermore, the determinant of the 
Poisson bracket between the gauge conditions and the gauge generators has a block 
diagonal form, so that it factorizes. Hence we obtain the integration measure as 

YlVNV7rNVN^V7riVV,,Vpr,VV^VpiVhijV7r^WXoVplVXVpx 

x5(p^)5(CA)5(7riv)5(vr.)5(a^)5(p0<^(^')<5(K)'^(b^)'5(^^')<5(PA) 

X v7I|det{x^K}l, (3.62) 


14 












where we denote = {'H'rp,'Hi). Integration over the variables N,7iisf, N\7ii,Vn,Pn, 
V'',Pi,X,px can be performed using the constraint ^-functions. The Hamiltonian and 
momentum constraints reduce to 03.541) and 03.551) . respectively, and we denote them 
collectively as = ("Ht, 'Hi). In addition, we write as an integral over auxiliary 

helds N'^ = {N,N^), essentially reintroducing the lapse and shift functions. The 
auxiliary helds are displaced so that the gauge hxing functions (/, /*) are cancelled, 
N + f ^ N and + /* —)■ iV*. Then the path integral is written as 


^dug=AA-i f l[VNVN^Vh,jV7r^WXoVpl6{xn\det{x^,n.}\ 


X exp 


fi/* 


+ p^dtXo - NUt - N^'Hi) - Hs, 


, (3.63) 


where A/" is a normalization factor. Integration over p\ gives a 5-function that imposes 
dtXo = 0. Therefore we decompose Aq to a constant component and an average free 
component Aq over time as 

Ao(^) = —JX + Ao(t), (3.64) 

K 

where J diX^ = 0. The integration over Aq is performed, which gives 


^dug=AA-i / \{VNVN^VhiPTi^^{x^) \det{x^,U,}\ 




X exp 


h 


dt { / {n^^dthij - NUt - N^'Hi) - Hi 


Bt 


'St 


where we have redehned 


Ut = + ^A. 

^/h K K 


, (3.65) 


(3.66) 


Unlike in [1^, the path integral fl3.65p does not include integration over the cosmolog¬ 
ical constant A, since we presume that the boundary conditions of the path integral 
dehne the (asymptotic) boundary values of all variables, including the boundary value 
of A. In particular, the value of A is set to a constant both on the initial Cauchy surface 
and on the spatial boundary. An extension of the path integral with an integral over 
A will be considered in Sect. 13.3.21 

Then we perform the integration over the momentum 7r*T Since the measure is at 
most linear in the momentum, the integration is Gaussian, and hence the integration 
can be performed in the standard way (see e.g. 30j). The integration amounts to 
expressing the momentum as 

Vh 


TT 






(3.67) 
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and including the factor iV 2 This result into 
^dug=AA-i [ l[VNVN^VhijN-^h-h{xn\det{x^M\ 




X exp 


i /1 


h \ K 


dt / NVh - 2A) + 


'St 


. (3.68) 


Then we express the field differentials as 

= 2hNVNVN^Vhij, (3.69) 

and write N~*h~^ = Ng^^{—g)~^, and obtain the path integral as 

^DUG =^f~^ [ YlRg^ii^g^^i-gy^Ndix^) |det{x'','Htt}|exp (^SEH[g^,u, A]) , 


xl^ 


(3.70) 

where S'Enfs'/ti/, A] is the Einstein-Hilbert (EH) action with an unspecihed cosmological 
constant A, 

SEn[giM^,A] = - [ {R-2A) + - (f d^x^/\^\K:. (3.71) 

Jm ^ JdM 

In summary, the difference compared to GR is that the value of the cosmological 
constant A is included in the initial and boundary conditions, rather than being a 
coupling constant of the Lagrangian. 

The next step is to express the gauge fixing factor of fl3.70p in a more useful 
form. For that purpose we consider specihc gauge conditions for the Hamiltonian and 
momentum constraints. The present theory has the advantage of enabling the use of 
the same gauges for the diffeomorphism symmetry as in GR. 


Dirac gauge First we consider the Dirac gauge [31 

Xb = hijTT^^ ^ 0, Xd = dj 


0 . 


(3.72) 


We dehne an operator in terms of the gauge transformation of the gauge condi¬ 
tions fl3.72p as 

QbuC = {a^,^ Rucj , = (e,e) • (3.73) 

Evaluating the Poisson brackets, we obtain the components of the operator as (up to 
the Hamiltonian constraint) 


= + 3A) , 

Kj 

Q^i ~ ~ ^iAd) 


= -2h3 ( - -R^K ) dj - 2d, 




1 \ 9 

A Xi uklB O I ^ hika Q \ xi ^ .k ^ i 


Qd, = - b^^xWk + xA'd^. + ^.aId, 


(3.74) 

(3.75) 

(3.76) 

(3.77) 


®The factor comes from 


det 2 , where the de Witt metric Gijki is 

regarded as a symmetric 6 by 6 matrix with indices (ij) and (kl) ranging over the six unique compo¬ 
nents. 
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where the momentum tt*-^ is written in terms of the metric fl3.67p and we denote = 
In order to obtain a gauge invariant form for the functional determinant 


32| . an extra factor N is included into the components Hence we replace fl3.74p 


with 


Qdo = + 3A) , (3.78) 

K 

and fl3.75p with = —Nx^^i — NdiX^- For practical applications, the components 
of the operator could be simplified by using the constraints, in particular the 
gauge conditions fl3.72p . and even further using the quasiclassical approximation (see 

i3)- 

Finally, the path integral can be written as 


^DUG 


= M- 


-1 


JJ Vg^^Vr]pVc*VcVc*Vdg^^{-g) 2 exp 


i 

h 


[g^iu, A] 


- / {gpXD + c*QdoC + + c-QdqC + c*Q[)jC>) 


IM 


, (3.79) 


where we have introduced pairs of anti-commuting fields c,c* and c*,c*, commonly 
referred to as Faddeev-Popov ghosts (and anti-ghosts), and an auxiliary field for 
each gauge condition Evidently, the full expression for the action is noncovariant 
in the Dirac gauge. 


Transverse harmonic gauge: 


0 ; = V^g 


in' 


(3.80) 


Transforming to this covariant gauge is achieved via the Faddeev-Popov trick in the 
same way as in GR. The operator corresponding to this gauge is again obtained from 
the gauge transformation of the gauge conditions 03.801) as 


= du - ropik • (s.si) 


Thanks to the gauge invariant form of the integration measure 
is obtained as 



, the path integral 


^DUG — A/" ^ 


2 exp 


i 

h 


A] 


(Tx 


-Vi^X 


^ + dpk 


{dpir'^ck 


r'^dpc^ 


rdpck) 


(3.82) 


JM / J 

where and c* are the Faddeev-Popov ghosts. 

Matter can be included similarly as in GR. For simplicity we assume that no 
extra gauge symmetries or constraints are involved. Finally, we define the generating 
functional by including external source and for the metric and the matter fields 
4/, respectively. 


^dug[J] = / X[Vgp,VgpVclVc^V^g^\-g)--^ 




h 




'M 


d^x i-VpX^ + dpCl {dpig^’^ck - ^^dpc’^ - 


. (3.83) 
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We have shown that the path integral for the fully diffeomorphism invariant uni- 
modular gravity fl2.13p has the same form as the path integral for GR with a cosmo¬ 
logical constant. The crucial difference from GR is that the value of the cosmological 
constant is set as a part of the boundary conditions for the path integral. 

The quantum effective action for the DUG theory can be dehned in the exact same 
way as for GR, since there are no extra conditions on the metric and the path integrals 
have the same form. 

3.3.2 Ng and van Dam form of the path integral 

The path integral for unimodular gravity can be extended by including an integration 
over the cosmological constant A. Then the path integral takes the following form: 

ZnvD = J d^{A)ZY)UQ{A), (3.84) 

where dfi{A) is an integration measure for A, and the path integral for unimodular 
gravity, Zdug(^); is given in f|3.82p with boundary conditions chosen to be consistent 



integral over Aq was assumed to include integration over both A and Aq due to the 


decomposition fl3.64p . Here we show that the path integral 03.841) follows straightfor¬ 
wardly from the canonical path integral of the action 02.13p . when we consider the 
vacuum state of the universe to be a superposition of the states corresponding to 
different values of A [^. 

We emphasize that 03.84|) is a quite different path integral compared to the one 
we derived above 03.82p . In the path integral 03.82p . the value of the cosmological 
constant is set as a part of the physical boundary conditions, which (together with a 
semiclassical matter distribution) dehne the vacuum state of the system. Including an 
additional integration over A means that we are integrating over different boundary 
conditions, i.e., vacuums. Below we attempt to justify the integration of A properly. 

Let I A) denote the vacuum state of the universe that is consistent with a given 
value of the cosmological constant A and with other relevant boundary conditions. The 
path integral that we have obtained for the fully diffeomorphism invariant unimodular 
gravity fl3.82p represents the vacuum transition amplitude 

(A|A) = Zdug(A). (3.85) 

We assume that transitions between vacuums are prohibited if the vacuums correspond 
to different values of A, i.e., the states |A) are assumed to be orthogonal, 

(A|A') = 0 if A ^ A'. (3.86) 

Furthermore we assume that the states are nondegenerate, i.e., there exists one state 
I A) for each value of A. The vacuum state of the universe is written as a superposition 
of the states corresponding to different values of A as 

\n) = f dAuj{A) |A). (3.87) 
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Choosing the weight function a; (A) dehnes which states |A) are included in the su¬ 
perposition. Now the vacuum transition amplitude is obtained in the form fl3.84p 
as 

= J dA|a;(A)|2(A|A) = J dp(A)ZDUG(A), (3.88) 

where the measure d^{A) is dehned by the weight function as 

d/i(A) = |u;(A)pdA. (3.89) 


Using the semiclassical approximation and then the stationary phase approxima¬ 
tion, it was argued in that the path integral fl3.84p for pure gravity is dominated 
by solutions whose cosmological constant A = 0. In the presence of matter 03.831) . the 
same argument was used in 1^ to obtain that the path integral 03.841) is dominated by 
the solutions of the Einstein equation whose cosmological constant is approximately 


A = 27iG -^^ (3.90) 

Jm^ 

where p is the energy density of a perfect fluid. This result was argued to imply that 
03.9np is the most likely value of the cosmological constant. It is intri guin g that using 
the present day energy density as an estimate for the average density one obtains 
a result that is surprisingly close to the observed value of A (the observed A/G being 
about three times the present average energy density). 

The result fl3.9Up is based on a hidden assumption that the given value of A is 
included in the vacuum state fl3.87p . It was assumed that all states |A) are weighted 
equally, |a;(A)p = constant. This corresponds to a total lack of physical boundary 
conditions regarding A, and then using the path integral for Ending the most likely 
value of A. This is an interesting argument, but speculative and conceptually prob¬ 
lematic. We indeed need information on the boundary conditions in order to estimate 
the average value of p over spacetime. Even if we accept fl3.90p as a valid estimate for 
the value of A in our universe, estimating the average of matter energy density over 
the whole spacetime is challenging, to say the least. 


4 Quantization of the unimodular gravity with a 
fixed metric determinant 

4.1 ADM decomposition of the action 

The gravitational part of the action fl2.6p is written in ADM form as 


^UG[iV,iV\h,,-,A,T] 



Nyjh 

K, 


- A (iV^/h - 6o) 

+ Sb + (4.1) 


A Hamiltonian formulation of an action of this form has been considered in ji^, and 
our following analysis is similar in several ways. 
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4.2 Hamiltonian analysis 

The momenta conjugate to N, iV*, and A are the primary constraints: 

ttat 0, TTi^ 0, px ~ 0. (4.2) 


The Hamiltonian is obtained as 

H= d^x {NUt + - eoX +vnt^n + v]>^-Ki +vxPx) + Hj3t, 


'St 


where the super-Hamiltonian is defined as 






kl 


\fh 


K 


+ \/hA, 


the super momentum is defined as 


= -2hijDkTi^\ 

vp^,v]^,vx are Lagrange multipliers, and Cq is the fixed scalar density. 

Preservation of the primary constraints implies the secondary constraints: 

■Hr ~ 0, Ri^O, U = N\/h — eo^O. 


(4.3) 


(4.4) 


(4.5) 


(4.6) 


The momentum constraint fl4.5p can again be extended with terms that are propor¬ 
tional to the primary constraints tin and px, 


'Hi = -2/ii,r>tir^'= + diNTTf, + dXpx, 


(4.7) 


since then it will generate spatial diffeomorphisms on S* for all the variables that are 
involved in the secondary constraints fl4.6p . The smeared Hamiltonian and momentum 
constraints fl3.28p satisfy the following Poisson brackets: 


{'Ht[^],'Ht[v]} = [ d^x{^dip - pdi^)R^{'Hj - O^Ntin - djXpx), 

mx%nT[^]} = nT[x%^], 

- ^P^d,x% (4.8) 

The Hamiltonian and momentum constraints have nonvanishing Poisson brackets with 
U: 

{U, WtKI) = (4.9) 

{W, <h[x*]} = x'^^iN^/h) + dix'NVh ^ eodixX (4.10) 

Hence the preservation of U, 

dtU = {W, H} ^ + eodiN^ + VhvN ^ 0, (4.11) 

2x/h 

is ensured by fixing the Lagrange multiplier as 

vm = wm = (4.12) 


20 








The preservation of "Hr, 


dtHT = {HT,H}^Vhvx^Q, 


fixes the Lagrange multiplier v\ as 
The preservation of Hi ~ 0, 


v\ = 0. 


dtUi = {H^, H} ^ eodiX ^ 0, 


(4.13) 

(4.14) 

(4.15) 


requires the introduction of the secondary constraint fl3.30l) . which was also present in 
the generally covariant formulation. The constraint fl3.30p is preserved in time since 
the Lagrange multiplier of the primary constraint pa ~ 0 has been fixed to zero fl4.14l) . 
We do not need any further constraints, but we still need to analyze and classify the 
existing constraints properly. 

We again decompose the variables A, pa as in fl3.4UI) and replace the constraint fl3.3Up 
with fl3.46p . The second class constraints A ~ 0 , pa ~ 0 can be used to eliminate 
the average free variables X,p\. Since the Hamiltonian constraint Ht contains the 
remaining zero mode Aq, and the zero mode p° of the primary constraint Pa ~ 0 
remains, we should also decompose Ht as 

'Ht= . + no= [ d^xHT, [ d^xnT = 0, (4.16) 

j^^dPxyh Jst JT.t 

where the zero mode and the average free component are respectively defined as 


\/h. 


Ho = d?x —-) + Aq / d^xVH ^ 0 (4.17) 

Jst \Vh 


and 

Ht 


K 

Vh 

K 

\/li 




Vh 

^ dPxVh JT.t Vh 


d^x ( 1 3. 


(4.18) 


Note that fl4.18p does not involve the cosmological variable Aq. The components Hq 
and Ht satisfy the Poisson brackets 


{•Ho,-Ho} = 0 , 

{Ho,Ht[^]}= f d^xdVRVRj-djNTTN-VXpV, 

J'£,t 

{HT[l],HT[ff\} = f d^xi^dip - {Hj - - djXpx), 

JEt 


where the smeared constraint Ht[^] is defined so that 


(4.19) 


HT[^]=mo + 'HT[^], Ht[^]= [ d^x^Hr, (4.20) 
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and where ^ is decomposed as any scalar, ^ and 

1 


eo = 


Ie, „ 'St 


(fxVh^, / d^xVh^ = 0. 


'St 


'St 


n,' n, , 

'Hr-HT + 


d^x {NH’t + N^H[) , (4.22) 


In the Hamiltonian fl4.3p . we obtain 

f d^x {NUt + N^Ui + wnt^n) = 

JEt 

where we have extended the Hamiltonian and momentnm constraints as 

0 , 

0 , 


n' = n^ + eo^t (^ 


and fnrthermore 


y/h 


d^xNn'T = No'H'o+ / d^xNn'T, 
JEt 


(4.21) 


(4.23) 


(4.24) 


where the zero mode and average free component of H'j, are dehned in the same way 
as for T-Lt in fl4.16p - fl4.18p . The hrst class constraints are the average free Hamiltonian 
constraint and the constraints "H', TTj. These constraints are associated with the 
invariance of the action fl2.6p nnder the metric determinant-preserving diffeomorphism 
fl2.ip - fl2.5p . Since the lapse N is not an nnspecihed mnltiplier in the Hamiltonian, we 
shonld add the term into the Hamiltonian density, where vt is an nnspecihed 

Lagrange mnltiplier. The Hamiltonian fl4.3p is rewritten as 


H= d^xiN'H'T-eoXo + VTn'T + Nm'. + v^Tii)+H b,. (4.25) 

JT:t 

The local constraints tt^v ~ 0, W ~ 0 and the zero mode constraints ~ 0, "Hq ~ 0 
are the second class constraints. The second class constraints tin ~ 0,W ~ 0 can be 
nsed to eliminate the variables N, tin as 

iV = ^, TT^ = 0. (4.26) 

The zero mode constraints ~ 0, "Hq ~ 0 can be nsed to eliminate the variables Ao,p° 
as 

The Dirac bracket that corresponds to the second class constraints {tin ~ 0,W ~ 
0,p° ~ 0,770 ~ 0) can be shown to be eqnal to Poisson bracket for all the remaining 
variables. 

When the second class constraints are set to zero strongly and the anxiliary vari¬ 
ables are eliminated as fl4.26l) and fl4.27l) . we obtain the Hamiltonian as 


H 


k f ^3^ 

d'^x\fh JEt 


K 

7^ 




kl 







+ HBf 


(4.28) 
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4.3 Path integral 

The canonical Hamiltonian for the action fl2.6p is written as 

H^= [ £x{NnT + N^'Hi-eo\)+HB^, 
Jst 


(4.29) 


where the Hamiltonian and momentnm constraints are defined in 04.41) and 04.51) . The 
second class constraints are 


U = NVh — eo ~ 0, ttn 

A 0, pxr 


Ur 


The hrst class constraints are vTj 0 and 


0 , Pi 


Hrp — PLt T 

n = n 


K tin 


2 

^odi 


\/Tl \/Tl 
Tv 
x/h 


iO, 

0 , 

0 . 


0 , 


(4.30) 


0 . 


(4.31) 


We denote the latter two constraints collectively as = {'H'rp,'H[). 

The gange fixing condition for W is defined as in 03.6Up . bnt there is no gange 
condition for N dne to the first pair of second class constraints in 04.301) . The gange 
conditions read as 

a* = W-r^0, (4.32) 

where one of the conditions has to be average free, so that the nnmber of gange 
conditions matches the nnmber of generators exactly. We choose it to be the zero- 
component, since the zero mode of the snper-Hamiltonian is a second class constraint, 
and hence we denote = (x°, x*)- 

The generator and the gange condition x° both snffer from a nonlocal linear 
dependence over the spatial hypersnrface, since their spatial inte gral s vanish by def¬ 
inition. The proper treatment of linearly dependent generators 28| is discnssed in 
Appendix 

The canonical integration measnre for the path integral is written as 


Y[VNV7iNVN^V7iiVhijV7i^WXoVplV\VpxS{U)6{7iN)S{\)S{px 










det 


(4.33) 


The initial and bonndary conditions on the cosmological variable are similar to DUG, 
i.e., the valne of A is set to a constant on the initial Canchy snrface and on the spatial 
bonndary. When integration over the variables N, ttn, T, A, px and pi is performed 
nsing the constraints, we obtain 




det 


X exp 


'St 


v/h' 


^ I dt{ I ( TT^^dth^j - -^Ut - n + epAp ) - Hs, 


, (4.34) 
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where we denote l-ii, = {TIt, 'Hi) and the Hamiltonian constraint is given in fl3.54p . Ex¬ 
pressing the ^-functions 5{Hq) and 5{H^) in terms of integrals over auxiliary variables 
N = (A^O) H) and and shifting the variables as N ^ N — ^ and A^* —>■ A^* — 
we obtain 


ZuG J Yl'I^NVN'^VhijVTi^WXo 


mn 




X exp 


n 


'St 


{ir^^dthij - NHt - N^Hi + 6 oAo) - 


Integration over the momentum vr*-^ gives 


(4.35) 


j]\Vg^„V\„g«\-g)-i (^N ^ 6(r) 


det \x'',H, 


X exp 


T ( SEulg/iu] — dt Ao (\/^ — Co) 


(4.36) 


where SEuidfn] is the Einstein-Hilbert action without a cosmological constant. Since 
the zero mode Aq depends only on time, integration over this variable gives a ^-function 
that imposes the unimodular condition fll.2p to hold on each slice of spacetime in 
average, 


^ug=AA-i f Y[Vg,,g^%-g)-h 

'' xt^ 




det <x^,Hg 


Hr 


X exp ( -S'eh[5' flV ] ] ■ (4-37) 


The integrated unimodular condition in the above path integral. 


[ {V^ — eo) — 0, 

JSt 


(4.38) 


does not constrain local deviations from the unimodular condition fll.2p as long as the 
average value of \/^ over S* remains fixed to that of eo. This is a quite surprising 
result, since we expected to see the unimodular condition to be imposed locally, like 
in the path integral for the HT action 1^. On the other hand, it makes some sense 
that quantum fluctuations around the classical held equation fll.2p are permitted. The 
physical purpose of the condition fl4.38p is to ensure that the number of physical degrees 
of freedom in the path integral fl4.37p matches that of DUG and GR, since together 
the gauge conditions r ^md the condition fl4.38p impose four conditions per point in 
space. 

In this theory, the quantum effective action is a function of the perturbative grav¬ 
itational held which satishes an integrated condition. Namely the trace of the 
perturbative held must have zero integral over at all times. 


n = 0 . 


(4.39) 


'St 


^Note that the change of variable N ^ N' = N + has a unit Jacobian, despite the fact that 
the transformation involves y/h. 
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In other words, the quantum effective action is built in the same way as in the HT 
theory [l^ except that the condition on the perturbative gravitational held (/^ = 0) 
is replaced with the integrated condition fl4.39p . The gravitational field is further 
constrained by the gauge conditions 

Counting of physical degrees of freedom 

In both cases, DUG and UG, Dirac’s counting of physical degrees of freedom gives the 
same result: two propagating modes plus one zero/single mode. In DUG, the extra 
zero mode is the cosmological variable A, which is a constant spatially and does not 
evolve. Thus the extra zero mode is not a true physical degrees of freedom. It is just a 
cosmological constant. Hence the physical degrees of freedom are the same as in GR. 

In the UG theory with fixed metric determinant, the canonical structure is partially 
different from DUG. In particular the integral of the Hamiltonian constraint is a second 
class constraint, and hence the gauge/coordinate conditions must contain one zero 
mode less than in DUG and GR. This may appear to imply that the extra zero mode 
would be a physical degree of freedom, but our analysis shows otherwise. This is 
evident in the path integral fl4.37p . where the single 5-function eliminates one zero 
mode by imposing the integrated unimodular condition fl4.38p . In other words it 
acts like an extra gauge/coordinate condition, so that the total number of conditions 
matches DUG and GR. Thus the actual number of physical degrees of freedom in UG 
is the same as in DUG and GR. 


Gauge fixing example We can choose the gauge conditions, for example, as 


Xu = (In h — $) = vTi l^ln h — d* — -j—Vh (In h — d*) J ~ 0, 

x\j = dj ^ 0, (4.40) 

where d) is a fixed function, and we denote the conditions collectively as Xu = (Xm Xu)- 
The first gauge condition Xu fixes the average free component of Ink. That is the 
average free component of the first condition of the Faddeev-Popov gauge j^. The 
gauge conditions xjj harmonic conditions on the spatial hypersurface. 

We define an operator Qu^, in terms of the gauge transformation of the gauge 
conditions (14.40^ as 




r = (f.e) ■ 


We obtain the components of the operator as 

Ql, = 2VhK + tjK, 

Qui = 2VhDi + Xu-Dj) 

Q\jq = -2^/h dj - 2dj Vh 

= -VhS)h^%di - 6)x’hdk + X\jd, + d,x\j- 


(4.41) 


(4.42) 

(4.43) 

(4.44) 

(4.45) 
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The components fl4.4211 - 04.4511 of the operator could be simplified by using the con¬ 
straints and in particular the gauge conditions. Finally, the path integral is written 
as 


[ Y[Vg^,Vr]VTi,Vc*VcVc*Vdg^\-g)-h 

•J 


exp 


h 


SehW - [ + + ?0?iciC + ?(3u.c* 

Jm ^ 


, (4.46) 


where we have introduced pairs of anti-commuting ghosts c, c* and c*, c*, and auxiliary 
helds fj^gi for each gauge condition. The helds c, c*, fj have vanishing average over 
space, since they are associated with the generator Kr and the gauge condition 
Including matter helds and dehning the generating functional can be done similarly as 
in 03.8311 . Evidently, the above expression for the path integral is not covariant. The 
presence of integration over space in both the averaged unimodular condition and the 
dehnition of average free helds renders the expression noncovariant. 

It indeed appears to be impossible to cast the path integral 04.3711 into a fully 
covariant form. The underlying reason is the fact that the zero mode of the super- 
Hamiltonian is a second class constraint, and hence one of the gauge conditions must be 
average free over space. In order to achieve a covariant description, we have to enlarge 
the gauge symmetry so that the total super-Hamiltonian becomes a gauge generator. 
This was achieved in Sect. |3l where a generally covariant form of unimodular gravity 
is considered. 


5 The canonical relation of the two theories 

In the case with a hxed metric determinant, it is crucial to notice that the Hamil¬ 
tonian fl4.3p is not a constraint, since it contains the term — d^xe^X. Therefore 
the bulk part of the Hamiltonian does not vanish on the constraint surface. This is 
a striking difference compared to Hamiltonian of the fully diffeomorphism-invariant 
theory 03.2011 . which is a sum of first class constraints. However, there exists a clear 
relation between these Hamiltonians, since the nonvanishing term can be eliminated 
(or introduced) via a simple time-dependent canonical transformation. 

Consider the following two canonical transformations of the variable p\ 

Pa = p'x ± eoC (5.1) 

with all other variables remaining unchanged. These two transformations are generated 
by the functionals 

F±= f d^x (Apa ± eoXt) , (5.2) 

respectively. The Hamiltonian transforms to 

H' = H+ ^ = H± d^xeoX. (5.3) 

dt ist 
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We can see that the transformation generated by F+ eliminates the nonvanishing term 
from the Hamiltonian fl4.3p . while the transformation generated by F_ introduces the 
nonvanishing term into the Hamiltonian fl3.2Up . Notice that the variable p\ appears 
only in the primary constraints Ca ~ 0 and pa ~ 0 of the two theories, and these 
constraints drop out of the Hamiltonian due to the consistency conditions for their 
Lagrange multipliers fl3.35p and fl4.14p . respectively. 

The theory with hxed metric determinant can be shown to be a (partially) gauge 
hxed version of the fully diffeomorphism-invariant theory. When we introduce the 
following gauge hxing conditions into the Hamiltonian fl3.48p FI 

hi = nVH — Co ~ 0, Pa ~ 0; Pa ~ 0, ^ 0, (5.4) 

and together with the second class constraints Ca ~ 0 and pn ~ 0, we obtain a Hamil¬ 
tonian that has the same form as 04.281) . except for the extra nonvanishing term in 
04.28p . — fj.^(Pxeo\o with Aq given in 04.27p . That extra term can be introduced 
with the canonical transformation 05. ip of the variable p\. Thus the theory 02.6p 
is a (partially) gauge hxed version of the theory 02.13p . In other words, the fully 
diffeomorphism-invariant theory dehned in 02.13p (and analyzed in Sect. [3]) is a gen¬ 
eralization of the unimodular theory of gravity with an enlarged gauge symmetry. 


6 Conclusions 


We have studied path integral quantization of two versions of unimodular gravity. 
In the fully diffeomorphism-invariant theory dehned by the action 02.13p . the path 
integral has the same form as the one of GR with a cosmological constant A 03.821) . 
except that the value of A is not set by the action. The cosmological constant A 
is an unspecihed value of the variable A. There exist two approaches regarding the 
interpretation of A in this theory: 


(i) The value of A can be set in the boundary conditions of the path integral, since 
it is a boundary value of the variable A. In this case, the value of A is completely 
unspecihed by the theory, and hence it needs to be set to the desired value 
by hand. One can use anthropic arguments for limiting the range of possible 
values of A (see jl-I^ for reviews), but we do not consider such arguments here. 
Physically, it makes no diherence whether the observed value of A is hxed by 
the boundary conditions or by setting the value of a coupling constant in the 
Lagrangian. Thus this approach is physically equivalent to GR. 


11 


Since the value of A is unspecihed, the vacuum state of the universe can be 
dehned as a superposition of vacuum states corresponding to diherent values of 
A [^. Such an approach was used in [l^, 1^ where the path integral of the 
form 03.841) was conjectured. A similar path integral was later obtained in 


14 


Starting from the action 02.131) . we have derived the path integral 03.841) without 
any addition or manipulation of variables. The integration over A arises due 
to the dehnition of the vacuum state 03.871) . The given theory shows that it is 
unnecessary to impose the unimodular condition on the metric determinant in 
order to obtain the path integral 03.84|) . 


®The first two gauge conditions are associated with the first class constraints ttn 
(the zero mode of Ht ~ 0), respectively. 


0 and Ho ~ 0 
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Using the semiclassical approximation and the stationary phase approximation 


one can argue 1^ that the path integral fl3.84l) is dominated by the values of 


A around the average energy density of matter over spacetime fl3.9Up . It is pre¬ 
sumed that the given values of A were included in the vacuum state fl3.87p . This 
result is interesting but problematic. In order to estimate the average energy 
density of matter over spacetime, we need information on both the matter and 
gravitational (background) helds, which depend on the assumed value of A. It 
could be interesting to search for alternative mechanisms that would single out 
the most likely values of A within the fully diffeomorphism-invariant theory. 


In the more conventional case dehned by the action fl2.6p . the path integral fl4.37p 
differs from the path integral of GR in two ways: (i) since the zero mode of the super- 
Hamiltonian fl4.17p is a second class constraint, the hrst class Hamiltonian constraint 
fl4.18p and an associated gauge condition have zero average over space, and (ii) the 
metric in the path integral must satisfy the integrated unimodular condition fl4.38p . 
The condition fl4.38p imposes the unimodular condition fll.2p to hold in average over 
space at each moment in time. The path integral has a generally noncovariant form 
due to the given differences. The perturbative gravitational held in the (semiclassical) 
quantum effective action must satisfy the integrated condition fl4.39p . 

At quantum level the unimodular condition can manifest itself in three ways. In 
the HT theory j^, the unimodular condition is imposed locally in the path integral and 


in the quantum effective action [l^. In the path integral and the quantum effective 


action of the UG theory fl2.6p . the unimodular condition is averaged over space fl4.38p . 
Lastly, the DUG theory fl2.13p does not involve a unimodular condition. 

In Sect. El we established the canonical relation of the two considered versions of 
unimodular gravity. While the actions fl2.6p and fl2.13p are shown to be equivalent 
classically, the time-dependent canonical transformation (15.1 p involved in the relation 
of their Hamiltonian structures has an interesting effect to the quantum theory. That 
is the appearance of the averaged unimodular condition (I4.38p in the path integral of 
UG (14.461) . Furthermore, the gauge symmetry is restricted, since the integral of the 
super-Hamiltonian over space (14.171) becomes a second class constraint. This implies 
that the path integral involves a pair of ghost helds and a Lagrange multiplier held 
whose average values over space must vanish. In practice, both of these implications are 
inconvenient to work with. Thus the fully diheomorphism-invariant theory considered 
in Sect. [3l or the previously worked out HT theory, are the preferable versions of 
unimodular gravity for quantization. 

The diherences in the path integrals of diherent versions of unimodular gravity do 
not necessarily imply that the physical predictions of the theories are diherent. The 
DUG and HT theories can indeed be expected to be physically equivalent, since the 
theories are related by a simple change of an auxiliary variable (see below (l2.13p L 
However, in addition to gauge hxing, the canonical relation between DUG and UG 
involves the time-dependent canonical transformation (15.ip . which leads to the afore¬ 
mentioned complications. Therefore it is still unclear whether the path integrals (I3.82p 
and (I4.46P produce equivalent predictions. Gonhrming this would require the formu¬ 
lation of Feynman rules and the calculation of the scattering matrices. This is a very 
demanding task in itself, which we wish to investigate in further work. 

Proper quantization of gravity requires more advanced methods. Two known ap¬ 
proaches are the spin foam models and the dynamical triangulations. Some steps 

















































toward loop quantization of unimodular gravity have already been taken in [26|], and 


more recently in [20 
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A Appendix: Quantization of gauge theories with 
linearly dependent generators 


The unimodular theories of gravity involve certain local gauge generators whose in¬ 
tegrals over the spatial hypersurface vanish by dehnition. This type of a spatially 
nonlocal linear dependence of generators is an inherent feature of the unimodular 
gravity theories, where the cosmological constant appears as a (constant) value of a 
scalar variable. Quantization of gauge theories with linearly dependent generators was 
achieved in j^. The Batalin-Vilkovisky formalism is suitable for the description 

of the nonlocally linearly dependent generators of unimodular gravity. In particular, 
the average free nature of the Faddeev-Popov ghosts and auxiliary fields associated 
with the average free generators is explained naturally within the given formalism. 

First we review the Batalin-Vilkovisky formalism. The formalism was applied to 
the minisuperspace formulation of Friedman-Robertson-Walker cosmology models in 
( 2 ^ , where a review of the formalism for theories with only bosonic gauge helds is also 
presented. Since the gravitational sectors of the unimodular gravity theories involve 
only bosonic helds, our presentation follows 29| with a few conventional differences due 
to the following application to the unimodular gravity theories considered in sections E] 
andm 

When the generators Ga are linearly dependent, there exist right zero eigenvectors 


Z a 
a 5 


GaZ: = 0 . 


(A.l) 


Here the condensed index a labels each local generator at every point in the spatial 
hypersurfaces. Hence summing over such an index involves an integration over space in 
addition to a sum over the components]^ The Latin index labels the zero eigenvectors 
a = 1,..., A. Here the vectors are linearly independent, i.e., we consider a hrst- 
stage reducible theory. The gauge conditions y" have to be similarly redundant as the 
generators, so that there exists left zero eigenvectors 



(A.2) 


The eigenvectors Z" and Z“ are the right and left zero vectors of the degenerate 
Faddeev-Popov operator, 

= (A.3) 

®Unlike in the sum over a condensed index does not involve integration over time. Furthermore 
we do not consider Euclidean quantum gravity, i.e., Wick rotation of time is not performed. 
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respectively. Thus, within this formalism, the Faddeev-Popov ghosts c", c* become 
gauge helds that require additional gauge hxing. For that purpose the set of Lagrange 
multipliers and ghosts (c", c*, T]a) is extended to [isj l 


$g = (c“,c;,r7„,c“,c':,i?“,0„r). 


(A.4) 


where c“, c*, 9a, '9°' are Grassmann anti-commuting variables and the rest are com¬ 
muting variables. The path integral and the corresponding effective gauge hxed action 
are written as 


Z = j P0'P7r,P<hgexp , 

SeE = S-I dt[clQy^ + + Vaix^ + 


where 0* and tt* are the gauge helds and their canonically conjugated momenta, and 
S is the action without gauge hxing. The extra Lagrange multipliers {da, 'd“) impose 
the gauge conditions a;“c" and c*(t“ on the Faddeev-Popov ghosts, where the gauge 
parameters (a;“, cr“) are arbitrary. The variables C* and are the ghosts for the 
Faddeev-Popov ghost helds. The so-called extra ghosts regulate divergent factors 
5(0) that appear in the original gauge hxing 5(x") with a redundant set of gauge 
conditions flA.2p . 

Integration over the ghost sector gives the path integral as 

^ (i^) . (^-6) 

where the gauge hxed Faddeev-Popov operator is dehned as 


X'cv _ /OQ: I ,a 

d' Q— 




(A.7) 


and the following matrices are introduced 


, A rya ya _ ya : 


(A.8) 


The path integral flA.OD is independent of the chosen gauge parameters (a;“, a"), since 
both the ratio of determinants (det detg^detg^) and the re gula ted gauge hxing 
factor are invariant under a change of the gauge parameters (see for a proof). 

Next we apply this formalism to the quantization of the two unimodular gravity 
theories (DUG and UG). 


A.l Fully diffeomorphism-invariant unimodular gravity 

Let us consider the quantization of DUG presented in Sect. [3l The second class 
constraints are given in fl3.58p . The generators are 

Ga = [TTN,'n'i,Pi,'HT,911,0'] (A.9) 

with fl3.59p . Gauge hxing conditions are chosen as in 03.601) . 

X" = ■ (A.IO) 
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The generator C and the corresponding gange condition p\ exhibit a nonlocal linear 
dependence, since their integrals over space vanish by definition. Hence there exist a 
single right zero vector, 


= 

= 


0 , 0 , 0 , 0 , 0 , 




\/h 


d?x\/h 


d?x\/liC' = 0, 


(A.ll) 


and a single left zero vector as 


Z„ = [0,0, 0,0, 0,1], 

r , _ (A.12) 

= / d^xpx = 0 . 

Jut 

Since only one pair of zero vectors exist, we have dropped the label a from the zero 
vectors and also from the other variables of the path integral flA.Op . 

We choose the gauge fixing parameters for the ghosts as 


= [ 0 , 0 , 0 , 0 , 0 ,- 1 ], 


a = 


0 , 0 , 0 , 0 , 0 , 




L (PxXh 


(A.13) 


Thus the ghost fields (c, c*) associated with the generator C are imposed to satisfy the 
gauge conditions 


CUaC" = — (TXC = 0, 
1 




d^xVhc* = 0 . 


We obtain the (now one-dimensional) matrices flA.811 as 


q = UaZ°‘ = — 1 , q = Za<J°' = 1 . 


(A.14) 


(A.15) 


In the amended Faddeev-Popov operator flA.711 . the gauge fixing term with flA.lOp 
contributes to the part of c* that involves the ghosts c, c* as 



d^x£yc*{x) {px{x),C'{y)] c{y) + = 



(A.16) 


This implies a unit contribution to the canonical measure of the path integral. 

The gauge condition p\ = 0 imposes p\ to become proportional to a spatial constant 
Pa (1230]), which is the integrated value of over space fl3.4ip . The extra ghost E 
introduces an independent term into this gauge condition, so that the integral of the 
condition over space no longer vanishes. 


(x“ + a^E) = E, 


(A.17) 
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which serves as a regulator for the corresponding ^-function in the path integral. The 
functional integral over E forces the Lagrange multiplier of the gauge condition px to 
have vanishing average value over space, 


VrjVE exp 


d'^xr] I Px + . ^ rr ^ 
IM \ d-^xy h 


oc 




Ie, ) 


’M 


d^xppx , (A. 18) 


which prevents the appearance of divergent factors 5(0). In sections and IT731 every 

5-function for an average free constraint is regulated in this way, and we denote such 
5-functions simply as 


Bijexp I ^ I d*xr)px') =d(px), 


(A.19) 


where the auxiliary field p is now assumed to have vanishing average over space. 

Once the additional gauge fixing flA.14D on the ghosts associated with the linearly 
dependent generators and gauge conditions is performed, and the 5-functions of average 
free constraints are regulated, it is easy see how the canonical path integral is obtained 
as fl3.63p after the nonphysical variables in the canonical integration measure fl3.62p 
have been integrated out (except for those variables deleted by the gauge conditions 
that are unspecified). In summary, the path integral obtained in Sect. 13.31 corresponds 
to the specific choice of the gauge parameters flA.131) . which are responsible for the 
additional gauge fixing required by the linearly dependent generators. 


A.2 Unimodular gravity with a fixed metric determinant 

Here we consider the quantization of UG presented in Sect. HI The generators are 

G„= (A.20) 

with fl4.3ip . Gauge conditions are chosen as in 04.321) . 

x“=k‘.y.x‘]. (A.21) 

We again have a nonlocal linear dependence, since the integrals of T-L'rp and over the 
spatial hypersurface vanish. A single pair of zero vectors is obtained as 


= [ 0 , 1 , 0 ], = / d^xn'T = 0 , 

JEt 

Za = [ 0 , 1 , 0 ], = f d^xx° = 0 . 


(A.22) 


The gauge fixing parameters Ua and a°‘ can be chosen so that the ghosts (c, c*) 
associated with the generator l-L'rp satisfy the condition of vanishing average value over 
space. The parameters are chosen as 


UJrv — 


(J = 


0 ,- ^ ,0 

k,d^x^fh 

0 ,- ^ ,0 


(A.23) 
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and the gauge conditions on the ghosts read as 


UnC = 


* OL 

c^a = 


d^x\/h 


d^x\/lic = 0, 


d^xVhc* = 0. 


(A.24) 


The determinants of the matrices flA.Sp have unit values 

q = 0 JaZ°‘ = 1, q = Za(J°‘ = 1. 


(A.25) 


The second class constraints (14.30^ contain a pair of average free constraints (A ~ 0, 
Pa ~ 0), which have to be treated in a similar way as a nonlocally linear dependent 
generator and a gauge condition. In the path integral, the contribution of these con¬ 
straints is just a unit factor to the Faddeev-Popov determinant, which is quite similar 
to the case of the constraints C and p\ in DUG. Hence we shall omit the analysis of 
these constraints here. 

Lastly, we explain how the path integral in Sect. 14.31 is obtained from the present 
formalism. As was discussed above, all the 5-functions for average free constraints 
has to be regulated in order to avoid divergent 5(0) factors. Integration over the 
extra ghost E and the additional Lagrange multipliers (6*, P) produces the 5-functions 
that impose the Lagrange multiplier of the gauge condition and the ghosts (c, c*) 
associated with the linearly dependent generators to become average free over the 
spatial hypersurface, 

5(p„cr")5((:nac")5(c*(T“), (A.26) 


where the conditions for the ghosts are flA.24p and the condition for the Lagrange 
multiplier p is 


= 


d'^xy/h Jst 


/ d^xVhrj = 0 . 

J'St 


(A.27) 


In Sect. 14.31 the 5-functions flA.26p are omitted in the path integral fl4.46p . since 
the helds x^, c, c* are assumed to satisfy the conditions flA.24p and flA.27p from the 
beginning. The Faddeev-Popov determinant in the path integral is dehned in terms 
of ghost helds that satisfy the conditions flA.24p . so that no zero modes are present in 
the ghost sector that generates the determinant. In particular, the nontrivial part of 
the determinant in the path integral fl4.34p (and thereafter) is written ai^ 


det = / 'DE Vd' exp 


h 


dtcl \x^,nf,\d 


(A.28) 


where the hrst components of the ghost helds = (c, c®) and c* = (c*, c*) are assumed 
to be average free over space. Finally, we can see that the path integral obtained in 
Sect. 14.31 corresponds to the specihc choice of the gauge parameters flA.23p . which are 
responsible for imposing the necessary conditions on the ghosts flA.24p and on the 
Lagrange multipliers (lA.27p . 


^^Note that here the sums over the repeated indices fi, v include integration over the spatial hyper¬ 
surface. 
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